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TRANSVERSE  PROPERTIES  OF  FIBROUS  COMPOSITES 
P.  E.  Chen*  and  J.  M.  Lin** 


ABSTRACT 

Tne  transverse  stiffness  and  strength  :f  unidirectional 
fiber- reinforced  composites  have  been  calculated  by  using  the 
finite-element  method  and  the  von  Mises-I-iencky  criterion. 

Both  the  square  and  hexagonal  arrays  have  been  considered  for 
the  fiber  configuration.  The  conditions  of  perfect  bonding 
ana  tocal  debonding  have  been  included  in  the  strength  calcu¬ 
lations  .  Experimental  work  has  also  been  conducted  on  boron- 
alumanum  and  stainless  steel-aluminum  composites.  The  trans¬ 
verse  properties  of  such  systems  have  been  measured  as  functions 
of  fiber  volume  content.  The  theoretical  results  are  compared 
with  the  experimental  data  of  our  own  for -the  metal-matrix 
composites  and  those  of  others  for  glass-epoxy  composites. 


research  Dept.,  Monsanto  Co.,  St.  Louis,  Mo. 

•’..a-ceriais  Research  Lab.,  Washington  University,  St.  Louis,  Mo. 


INTRODUCTION 


Various  methods  [1-7]  have  been  proposed  for  the  calcu¬ 
lation  of  transverse  stiffness  of  fibrous  composites.  This 
paper  employs  the  finite-element  method  which  is  believed  to 
be  relatively  more  accurate  and  rigorous.  Perhaps  more  sig¬ 
nificantly,  the  present  paper  deals  also  with  the  transverse 
strength  of  fibrous  composites  using  a  method  based  on  the 
von  Mises-Hencky  criterion.  It  is  assumed  that  the  distortional 
energy  condition  is  valid.  This  assumption  is  substantiated 
by  the  fact  that  the  fibrous  composites  usually  fail  at  a 
relatively  low  level  of  transverse  loading,  and  also  by  the 
actual  stress-strain  behavior  of  the  materials.  Two  limiting 
conditions  have  been  considered  in  the  strength  calculations. 

Por  the  first  condition  the  fibers  are  assumed  to  be  perfectly 
bonded  to  the  matrix,  while  for  the  second  condition  the  fibers 
are  assumed  to  be  totally  debonded  from  the  matrix,  thus  pro¬ 
viding  the  upper  and  lower  bounds.  It  is  also  assumed  that 
the  fibers  are  circular  in  cross  section  and  unidirectionally 
aligned.  Moreover,  the  transverse  properties  have  been  calcu¬ 
lated  for  both  the  square  and  hexagonal  arrays  as  the  fiber 
configurations  in  the  normal  plane.  The  theoretical  results 
are  confirmed  experimentally. 
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TKEORY 


The  stresses  and  displacements  of  the  composite  in  the 
elastic  region  are  calculated  by  using  the  finite-element 
method  [8-10].  This  method  utilizes  the  direct  stiffness  con¬ 
cept  which  considers  the  composite  system  as  an  assemblage  of 
idealized  elastic  elements  assumed  to  be  joined  together  at 
discrete  nodes.  By  adding  together  at  each  node  the  stiff¬ 
ness  coefficients  of  adjacent  elements  a  stiffness  matrix  for 
the  system  is  obtained.  This  stiffness  matrix  relates  the 
external  forces  acting  on  the  nodes  to  the  displacements  of 
the  nodes.  By  inverting  the  stiffness  matrix  one  obtains  an 
influence  matrix  which  gives  the  nodal  displacements  as  a 
function  of  the  external  forces  or  loads  acting  on  the  system. 
Using  the  same  strain  pattern  for  the  elastic  element  that 
was  assumed  for  deriving  its  stiffness  coefficients,  one  may 
derive  a  matrix  of  stress  coefficients  which  gives  the  stresses 
in  the  element  as  a  function  of  its  nodal  displacements.  The 
method  is  ideally  suited  for  analyzing  multi-phase  materials. 

It  is  founded  on  equilibrium  and  compatibility  conditions. 

The  Typical  Region 

In  order  to  render  the  problem  more  tractable,  the  fibers 
are  assumed  to  be  ideally  packed  into  square  or  hexagonal  array 


as  shown  in  Figure  la  or  2a.  Further,  by  invoking  symmetry 
and  compatibility  conditions,  it  is  only  necessary  to  consider 
a  typical  region  as  shown  in  Figure  lb  or  2b.  The  procedure 
for  analyzing  the  typical  region  is  basically  similar  to  that 
described  in  Reference  11.  The  essential  details  in  the 
analysis  of  the  typical  region  are  given  in  the  Appendix. 

The  vcr.  "ises-Hsr.ckv  Crltericr. 

Various  theories  [12-15]  have  been  proposed  for  predicting 
the  strength  of  materials.  However,  the  theory  introduced  by 
von  Mises  [16]  and  reinterpreted  by  I-Iencky  [17]  is  generally 
recognized  as  conceptually  most  consistent,  and  it  is  also 
supported  by  experimental  evidences. 

The  total  strain  energy  stored  in  an  elastic  body  can 
be  divided  into  two  parts,  the  dilatational  energy  and  the 
distortional  energy.  The  von  Mises-Hencky  theory  postulates 
that  yielding  sets  in  when  the  distortional  energy  reaches  a 
critical  value.  For  a  uniaxial  and  plane  state  of  stress  the 
criterion  becomes 


a,2  -  a,a2  +  o22  =  sm2  (1) 

where  aj  and  are  the  principal  stresses,  and  Sm  is  the 
strength  of  the  matrix  material.  The  quantity  on  the  left 
side  of  equation  (1)  will  be  referred  to  as  the  normalized 
distortional  energy  hereafter. 
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The  Stiffness  and  Strength  Calculations 

To  proceed  with  the  calculation  of  the  transverse  prop¬ 
erties,  the  typical  region  is  first  divided  into  a  sufficient 
number  of  finite  elements,  the  procedure  given  in  the  Appendix 
is  then  used  to  calculate  the  distributions  of  stress  and  dis¬ 
placement  in  the  region.  Corresponding  to  the  assumed  displace¬ 
ment  conditions,  the  applied  normal  stress  can  be  calculated 
from  the  stress  distribution  previously  determined.  The  stiff¬ 
ness  is  obtained  by  a  simple  application  of  the  Hooke's  law. 

Assuming  perfect  bonding  or  total  debonding,  the  normal¬ 
ized  distortional  energy  is  evaluated  for  every  element  in  the 
domain,  thus  determining  also  the  maximum  normalized  distor¬ 
tional  energy.  For  the  condition  of  perfect  bonding  the  fibers 
are  assumed  to  be  in  perfect  contact  with  the  matrix,  and  the 
composite  is  considered  continuous,  from  the  mechanistic  point 
of  view,  at  the  interface  between  the  constituents.  For  the 
condition  of  total  debonding  the  fibers  are  assumed  to  be 
completely  separated  from  the  matrix.  The  transverse  strength 
of  the  composite  can  be  calculated  from  the  following  equation: 

st  =  sm  ac  /  (umax^  ^ 

where  a  is  the  applied  stress  on  the  composite  under  the 
assumed  displacement  conditions  as  described  in  the  Appendix, 
and  UJTlax  is  the  corresponding  maximum  normalized  distortional 
energy.  Both  ac  and  Umax  are  functions  of  fiber  volume  con¬ 
tent,  fiber  geometry,  condition  of  bonding,  as  well  as  the 
constituent  properties. 


mrr  - 
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EXP JRIMENT 


Materials 

Two  types  of  metal  have  been  used  as  the  matrix  materials , 
namely  the  6061  and  2024  aluminum  alloys.  Boron  and  stainless 
sceel  fibers  have  been  used  for  the  reinforcing  phase.  The 
mecnanical  properties  of  the  constituent  materials  relevant 
to  the  theoretical  calculations  are  given  in  Figures  3-6 

Technicue 


Solid  state  diffusion  bonding  technique  was  utilized  to 
fabricate  the  composites  used  in  the  experimental  work. 

The  general  bonding  conditions  were  3  hours  at  11,000  psi  and 
900°F  for  the  boron-aluminum  system,  and  30  minutes  at  6,000 
psi  and  900°F  for  the  stainless  steel -aluminum  system.  The 
equipment  used  for  diffusion  bonding  is  shown  in  Figure  10. 
Tensile  specimens  were  prepared  by  first  shearing  or  wet 
grinding  to  approximate  size  and  then  polishing  with  #180  SiC 
paper  under  running  water  to  final  size  of  0.30  in.  wide  by 
2.75  in.  long,  with  the  thickness  varying  from  0.030  in.  to 
0.116  in..  A  typical  tensile  specimen  is  shown  in  Figure  11. 
All  specimens  were  tested  by  an  Instron  machine  at  room  temp¬ 


erature  . 
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RESULTS  AND  DISCUSSION 


The  theoretical  approach  as  described  previously  was  used 
to  calculate  the  transverse  properties  of  boron-6061  aluminum, 
stainless  steel-2024  aluminum,  E  glass-epoxy  and  S  glass-epoxy 
composites,  as  functions  of  the  fiber  volume  content.  The 
theoretical  results  are  compared  with  our  own  experimental 
data  obtained  for  the  metal-matrix  composites,  and  the  experi¬ 
mental  results  obtained  by  other  investigators  [3,18,19]  for 
glass-epoxy  composites,  as  shown  in  Figures  3-9. 

The  composite  transverse  properties  have  been  calculated 
for  both  the  square  and  hexagonal  arrays  as  the  idealized  fiber 
packings  in  the  matrix.  It  is  interesting  to  note  from  the 
calculated  results  that  the  square  array  gives  relatively  higher 
transverse  stiffness,  but  lower  transverse  strength  than  those 
of  the  hexagonal  array.  However,  the  stiffness  seems  to  be 
less  sensitive  than  the  strength  to  the  change  in  fiber  geometry. 
Strictly  speaking,  the  fiber  geometry  is  neither  square  nor 
hexagonal,  but  it  appears  to  be  closer  to  the  square  array. 
Perhaps,  it  should  be  mentioned  in  passing  that  all  curves  for 
the  square  array  terminate  at  Vf  =  78.5%,  and  those  for  the 
hexagonal  array  terminate  at  =  90.6%.  The  above-mentioned 
fiber  volume  contents  are  the  geometric  limitations  for  the 
respective  arrays. 

In  addition  to  the  fiber  geometry,  the  extent  of  debond¬ 
ing  between  the  matrix  and  the  fibers  also  has  significant 
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efrect  on  the  transverse  strength  of  the  composites,  as  can 
be  seen  in  Figures  4,  6  and  9.  In  order  to  see  the  actual 
conditions  of  bonding  and  debonding,  the  specimens  were 
studied  under  the  optical  microscope.  Three  pictures  taken 
ur^er  the  microscope  are  included  in  this  paper  for  illustra¬ 
tion.  Corresponding  to  the  experimental  data  as  shown  in 
Figure  4,  the  condition  of  bonding  of  boron  fibers  in  6061 
aluminum  matrix  is  shown  in  Figures  12  and  13.  Related  to 
the  experimental  results  given  in  Figure  6,  the  condition  of 
nmcnomg  or  scam-ess  scee_  in  2324  alum.,  nun  matrix  is  shewn 


in  Figure  14. 

It  is  also  important  to  point  out  from  the  results 
o_.i_i.ineo  tnus  rar  that  m  no  case  has  the  composite  transvers 
strength  exceeded  the  matrix  strength.  The  decrease  in  trans 


verse  strength  for  higher  fiber  content  is  basically  caused 
by  the  corresponding  increase  in  stress  concentration  [20]  . 
The  increased  transverse  strength  for  hexagon:  1  array,  as 
compared  with  the  square  array,  is  believed  due  to  a  more 
efficient  shear  transfer  device  similar  to  that  as  suggested 
by  Sadowsky  [21] . 
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NQMENCLATURE 


Strength  of  the  matrix  material. 

Composite  transverse  strength. 

Young's  modulus  of  the  fiber  material. 

Young's  modulus  of  the  matrix  material. 

Composite  transverse  stiffness. 

Poisson's  ratio  of  the  fiber  material. 

Poisson's  ratio  of  the  matrix  material. 

Maximum  normalized  distortional  energy. 

Fiber  volume  content. 

Rectangular  coordinates . 

Displacements  in  x  and  y-directions . 

Displacements  in  x  and  y-directions  for  Case  1 
in  the  Appendix. 

Displacements  in  x  and  y-directions  for  Case  2 
in  the  Appendix. 

Principal  stresses. 

Average  normal  stresses  in  x  and  y-directions 
for  Case  1  in  the  Appendix. 

Average  normal  stresses  in  x  and  y-directions 
for  Case  2  in  the  Appendix. 

Average  normal  stresses  in  x  and  y-directions 
for  Case  3  in  the  Appendix. 

Applied  normal  stress  on  the  composite  for  Case  3 

in  the  Appendix  (a  =  a",  under  the  assumed  condi- 

c  x 

tions) . 


Shearing  stress  in  xy-plane  parallel  to 
y-axis . 

Width  of  the  typical  region. 


or 
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APPENDIX 


The  determination  of  the  stress  and  displacement  distri¬ 
butions  in  a  composite  as  shown  in  Figure  la  or  2a  can  be 
accomplished  by  analyzing  a  typical  region,  as  shown  in 
Figure  lb  or  2b.  The  region  is  so  chosen  that  under  a  normal 
stress  at  infinity,  the  rectangular  region  after  deformation 
remains  rectangular.  The  finite-element  technique  and  the 
method  of  superposition  are  used  to  solve  the  problem  in  the 
following  steps,  assuming  that  the  applied  normal  stress  is 
in  the  x-direction: 

1.  Solve  Case  1  which  is  defined  by  the  following  bound¬ 
ary  conditions: 

t „  =  0  along  the  entire  boundary, 

xy 

u  =0  along  AO  (points  remain  on  the  y-axis  because 
of  symmetry) , 

u  =1  along  BC  (arbitrarily  specified  unit  displace¬ 
ment)  , 

v  =0  along  OC  (points  remain  on  the  x-axis  because 
of  symmetry) , 

v  =0  along  AB  (specified  displacement  condition) . 


The  displacement  field  thus  calculated  is  (u^Vj),  and 
the  average  normal  stresses  in  the  x  and  y-directions  are  a 


and  a 

yi 


respectively. 
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2.  Solve  Case  2  which  is  defined  by  the  following  bound¬ 
ary  conditions: 

xXy  =  0  along  the  entire  boundary, 

u  =0  along  AO, 

u  =0  along  BC, 

v  =0  along  OC, 

v  =1  along  AB. 

Che  displacement  field  thus  calculated  is  (u0,v2)  and 
the  average  normal  stresses  in  the  x  and  y-directions  are  a 
and  cv2  respectively. 

3.  Solve  Case  3  which  is  characterized  by  o’  =  0 ,  solu¬ 
tion  of  Case  2  is  multiplied  by  (-a  /  o^2)  and  summed  with 
that  of  Case  1.  Thus  the  corresponding  applied  normal  stress 
on  the  composite  is 

a. 


X: 


'VI  - 

aC  -  ax  *  axi  -  If-  °x2 

ay2 


(3) 


Likewise  for  the  stress  and  displacement  components. 
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